Third harmonic generation in graphene and few-layer graphite films 
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We observe optical third harmonic generation from graphene and few-layer graphite flakes pro- 
duced by exfoliation. The emission scales with the cube of the intensity of the incident near-infrared 
femtosecond pulses and has a wavelength that is one-third of the incident wavelength, both consis- 
tent with third harmonic generation. We extract an effective third-order susceptibility for graphene 
that is on the order of 10~^^ m'^/V^, which is comparable to that for materials that are resonantly 
excited, but larger than for materials that are transparent at the fundamental and third harmonic 
wavelengths. By measuring a set of flakes with different numbers of atomic layers, we find that the 
emission scales with the square of the number of atomic layers, which suggests that the susceptibility 
of graphene is independent of layer number, at least for a few layers. 



Graphene is a monolayer of carbon atoms arranged in a 
hexagonal two-dimensional lattice. It has a linear disper- 
sion relationship between energy, E, and wavenumber, k, 
E{k) = zLhvpk, where the Fermi velocity vp « 10^ m/s 
[see Fig. 1(b)] [iHl]- Since the material became readily 
available less than a decade ago it has been the sub- 
ject of extensive investigations [e.g., see Refs. 0and0and 
references therein]. Graphene exhibits a number of un- 
usual and remarkable transport properties that make it 
attractive for nano-electronic applications 0-01, includ- 
ing high mobilities j^. Isl. lol- Jl^ and nearly-ballistic trans- 
port at room temperature [5|,|6[; however, it is the optical 
properties of graphene that are of primary interest here. 

The linear absorbance of graphene is flat and approx- 
imately 2.3% across the entire visible spectrum [llf. IT^. 
Thus, graphene can be considered to be both highly 
absorbing and/or highly transmitting, depending upon 
one's point of view or application. Since a negligible 
fraction (< 0.1%) is reflected [111, 97.7% of the incident 
light is transmitted. In this sense the sample is certainly 
transparent. On the other hand, if one were to assign 
an effective absorption coefficient, a, to a monolayer of 
thickness 0.33 nm (even though it may be questionable 
to use macroscopic parameters, such as a, for such thin 
samples), it would be very large (a « 7 x 10^ /cm). As 
suggested by this large effective absorption coefficient, 
graphene interacts strongly with light. The strong broad- 
band nature of the interaction of light with graphene 
is consistent with its linear bandstructurc, where inter- 
band transitions of roughly equal strength arc available 
throughout the visible. The combination of broadband 
transparency and the high mobilities mentioned earlier 
makes graphene a promising candidate for use as a trans- 
parent conductor Il3l- la j in photovoltaic devices 171 Il8| 
and touch screens [1J|. The high broadband absorption 
{i.e., optical conductivity) and high carrier mobility sug- 
gest that graphene may find applications in a range of 
optically-controlled transport devices, such as broadband 



and ultrafast photodetectors optical modulators 

26 , |27[ , plasmonic devices □roadband polar 



irizers 



34| , and magneto-optical devices [35 



The strong broadband coupling of light with graphene 
also suggests that it might be an attractive nonlinear ma- 
terial for large bandwidth photonic applications. How- 
ever, second order nonlinearities in graphene are ex- 
pected to be weak because of symmetry considerations. 
For example, graphene mounted on a substrate has C^y 
symmetry, and therefore, the only allowed sources of sec- 
ond order responses are the interfaces [3^. By compari- 
son, for multilayer graphite, the bulk can contribute, de- 
pending upon the stacking arrangement (STj . Indeed, sec- 
ond harmonic generation has been investigated both the- 
oretically 3^- io] and experimentally 37 , 41 1 in graphene 
and few-layer graphite films. Graphene flakes mounted 
on oxidized Si substrates have produced second harmonic 
signals, but the signal from the graphene is weaker than 
that from the substrate [13, EH . For bilayers and few- 
layer graphite samples similarly mounted, the second har- 
monic signal is dominated by the graphite films, but it is 
weak [la Ell- 
in addition, the third order nonlinear optical responses 
associated with saturable absorption of graphene 42- 
44 1, few-layer graphite films 13, 3-13], and suspensions 
of graphene flakes [s^] have been measured and time- 
resolved using a variety of pump-probe configurations. 
Typically, in these experiments, the interband absorp- 
tion of the pump pulse produces noncquilibrium elec- 
tron and hole populations, which subsequently relax by 
carrier-carrier scattering, phonon emission, diffusion and 
recombination. The time-delayed probe pulse is used to 
measure the dynamic change in the absorption coefficient 
caused by Pauli blocking of the final and initial states by 
the electron and hole distribution functions, respectively. 
Consequently, the emphasis in these studies has been the 
determination of the carrier, phonon and recombination 
dynamics. Because the measured nonlinearities in such 
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experiments are dynamic, cumulative, and incoherent, (a) 
they are not well suited to defining and quantifying a 
third-order susceptibility per se, and little attention has 
been given to quantifying the magnitude of these nonlin- 
earities. 

By comparison, measurements of coherent third order 
nonlincarities are relatively rare 58-63l- Two-photon ab- 
sorption has been measured in bilayer graphite flakes [58| , 
and a third order nonlinear susceptibility with magni- 
tude, Ix^'^^l, on the order of 10~^^ m^/V^ was extracted, 
which is consistent with the quantum perturbation the- 
ory used by the authors. By comparison, in the same 
study [53, no two-photon absorption was observed in 
monolayer graphene, also consistent with the authors 



theory which predicted a \x 
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m^/V^. In con- 



trast, however, ballistic photocurrcnts have been injected 
into multilayer epitaxial graphene samples by quantum 
interference between two-photon absorption and one- 
photon absorption, which suggests that the strength of 
the two-photon absorption is significant in graphene (sij . 



Moreover, 
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s as large as 10^^^ m^/V^ have been 
reported for coherent four-wave mixing experiments in 
graphene [HO], and a study on coherent self-phase mod- 
ulation in dispersion of graphene sheets containing effec- 
tively about 100 layers deduced a monolayer Ix^'^-'l of the 



same order of magnitude [61 



Here we report observation of third harmonic gen- 
eration (THG) from graphene and few-layer graphite. 
The samples studied are graphene and few-layer graphite 
flakes produced by mechanical exfoliation and mounted 
on a 300-nm silicon dioxide layer grown on a 500-fim- 
thick silicon substrate. Flake thicknesses are deter- 
mined by their optical contrasts (Hij and atomic force 
microscopy. Third harmonic emission is measured us- 
ing the experimental geometry schematically sketched in 
Fig. [TJa). The fundamental near-infrared pulses (w, red) 
with an average power in the range of 1 - 10 mW are 
obtained from an optical parametric oscillator pumped 
by a Ti:sapphire laser at a repetition rate of 80 MHz. 
Each pulse is tightly focused onto the sample by a mi- 
croscope objective, and the third harmonic generated by 
the sample (3w, green) is collected in the reflected direc- 
tion by the same lens. Color filters are used to block the 
refiected w beam. An adjustable mirror (M) is used to 
direct the 3a; beam cither to an imaging charge-coupled 
device (CCD) camera, a spectrometer, or a silicon pho- 
todiode connected to a lock-in amplifier. The camera is 
used to rapidly locate a flake having the desired thickness 
and to spatially image the 3a; spots. The spectrometer, 
which is equipped with a thermoelectric cooled silicon 
CCD camera, allows the detection of weak 3a; spectra. 
The combination of the photodiodc, which has an inter- 
nal amplifier, and the lock-in amplifier provides for the 
continuous monitoring of the 3a; power as experimental 
parameters are varied. 

The 3a; signal produced by a graphene flake when it 
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FIG. 1. (a) Experimental scliematic: A fundamental (lu, red) 
pulse is normally incident on the sample. The third harmonic 
(3a;, green) is detected in the reflected direction by a CCD 
camera, a spectrometer or a photodiode connected to a lock- 
in amplifier, (b) Schematic diagram showing the linear dis- 
persion of graphene, indicating the one-photon absorptions of 
bj and 3a; and illustrating the energy level description of third 
harmonic generation. 



is irradiated by an to beam with horizontal linear polar- 
ization is summarized in Fig. [2] The spectrum of the 
emitted 3a; light produced by an uj beam with an av- 
erage power of 10 mW is shown in the inset of Fig. [5] 
The central wavelength of 575.5 nm is within 0.4% of the 
separately measured uj wavelength (1720.4 nm) divided 
by 3 {i.e., 573.5 nm), which confirms that we are mea- 
suring the third harmonic [see Fig. [IJb) for the energy 
diagram of THG] . When the laser spot is moved from the 
flake to the bare substrate, no 3a; signal is observed un- 
der the same conditions. Also, the second-harmonic sig- 
nals from the graphene flakes, the Si/Si02 substrate, and 
their interfaces, which were previously observed in sim- 
ilarly mounted samples (stI . |41| , are much smaller than 
the third harmonic and are not measurable under our 
conditions. In addition, by using an analyzing polarizer 
in front of the spectrometer, we find that the emitted 3a; 
radiation is linearly polarized along the same direction 
as the UJ pulse. 

The main panel of Fig. [2] shows the 3a; average power 
as a function of the lu average power. For convenience, 
the corresponding lu peak on-axis irradiance and pulse 
fluence are also shown on the top axis. The average 3a; 
power is calibrated by comparing the measured 3a; signal 
on the spectrometer to that produced by a laser pulse of 
known power at 560 nm incident on the substrate. In this 
way, the reflection of the substrate, the loss of all optics, 
and the sensitivity of the spectrometer are included in 
the detector calibration. The measured signal is consis- 
tent with a cubic dependence on incident uj power, as 
indicated by the comparison with solid line in the main 
panel of Fig. [5] and as expected for third harmonic signal. 
The small deviation from cubic behavior at extremely low 
power levels can be attributed to unremoved noise back- 
ground in the signal. 

In an attempt to quantify the strength of the third 
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FIG. 2. The average power of the third harmonic as a 
function of the average power of the incident fundamental 
(solid squares), and the spectrum of the third harmonic (solid 
squares, inset). The solid line is a cubic fit (x^^"* = 4 x 10~^^ 
m^/V^) of Eq. (1) to the data using parameters appropriate 
for our experimental conditions, as discussed in the text. 



harmonic signal produced by graphcne, we now extract 
an cfFectivc third-order susceptibility, x'-'^^ . by comparing 
the 3a; signal shown in Fig. [5] with the well-known re- 
sult for a thin bulk sample in the limit of an undepleted 
fundamental 1361: 



1 i3uj 
4 2n3„c 



(1) 



where c is the speed of light in a vacuum, d is the sample 
thickness, n^^^ is the index of refraction of the material 
at 3a;, and S^i^ (E^) is the magnitude of the 3a; (w) field 
and where each field is related to its irradiance by Ii = 
riieocEiE* /2. Furthermore, if the irradiance associated 
with each fundamental pulse is taken to be Gaussian in 
space and time. 
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with a diameter W and a pulsewidth r (both in full width 
at half maxima), the peak on-axis irradiance, /q, is re- 
lated to the average power, P^, by 



""^V4h^41n2 



(3) 



where / is the repetition rate of the laser. While the 
use of bulk concepts, such as susceptibilities and indices 
of refraction, to describe a single atomic layer is ques- 
tionable, this procedure allows us to assign an effective 
material parameter that allows comparison of THG in 
graphene to that in other well-known bulk materials. Us- 
ing values appropriate to our experiment [i.e., d = 0.33 



nm, w rig^ = 2.4, W = 3.5 /Ltm, r = 320 fs, A =1720.4 
nm, and / = 81 MHz], these expressions are used to ex- 
tract a value of 4 x 10"" mVV^ from Fig. ^ 
We note that the sample used for these measurements 
is shaped like a rectangle with a width ~ 1 /j,m and a 
length much larger than the focused laser spot diameter. 
Consequently, the overlap between fundamental focused 
laser spot and the sample is imperfect. The fit shown 
in Fig. [2] is corrected for this incomplete overlap. The 
fit shown in Fig. [2] is not corrected for multiple reflec- 
tions associated with the layered oxide-Si structure. A 
calculation based on light propagation through our spe- 
cific layered structure. j4ll |63| which incorporates third 
harmonic polarizations, indicates that the x^^^ value ex- 
tracted in this case is changed only by roughly a factor 
of two from that reported above, which is within our ex- 
perimental uncertainty. Therefore, for pedagogical sim- 
plification, we ignore those complications here. We also 
note that the upper limit on the fiuence (and irradiance) 
shown in Fig. [2] is determined by damage to the sample. 

We also investigate the nonlinear optical response as 
the sample evolves from monolayer graphene to multi- 
layer graphite by measuring the third harmonic from 
flakes of thickesses ranging from 1 to 6 atomic layers 
mounted on a single Si/Si02 substrate for an average 
power of 1 mW, as shown in Fig. |3l The data point 
labeled as having zero atomic layers represents the re- 
sponse of the bare substrate, and it confirms that any 
third harmonic generated by the substrate is negligible 
on our measurement scales. The solid curve in Fig. [3] is a 
fit of Eq. 1 to the data for an effective third-order suscep- 
tibility of Ix^^-* I 1 X 10"^^ m^/V^, and it has a quadratic 
dependence on sample thickness. Notice that this value 
is ^ 2.5 times larger than that extracted from Fig. [2l To 
the degree possible, the experimental conditions are held 
constant for the two sets of measurements; however, they 
are performed on different days, on different flakes, of dif- 
ferent sizes, mounted on different substrates. Given the 
cubic dependence of the third order interaction on the 
field and given the inherent uncertainties and variation 
in the experimental parameters, this represents reason- 
able agreement between the two measurements. For these 
reasons, |x 

I ~ 10-1^ mVV^ should be regarded as an 
order of magnitude estimate of the susceptibility. This 
magnitude is one order of magnitude smaller than those 
estimated from recent four-wave mixing [Hoj and self- 



phase modulation |61l | experiments, but more than two 
orders of magnitude larger than that from two-photon 
absorption experiment and quantum perturbation calcu- 
lation [HI]. Further efforts, from both theoretical and 
experimental sides, are needed to fully understand third 
order response of graphene. 

The deduce effective x^'^^ is several orders of magnitude 
larger than that for transparent materials [3^. This is 
because third order response is dramatically enhanced by 
the resonant nature of the intermediate and final transi- 
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FIG. 3. The average power of the third harmonic as a function 
of the number of atomic layers for an average fundamental 
power of 1 mW (solid squares). The solid line is a quadratic 
fit (x'^' = 1 X 10"^" mVV^) of Eq. (1) to the data using 
parameters appropriate for our experimental conditions, as 
discussed in the text. 

tions needed for the THG [s^ . The nonlinear response of 
graphene is also broadband, owing to its bandstructure. 
On the other hand, however, x^^-* of graphene is compara- 
ble to those associated with resonant transitions in other 
materials {e.g. metal nanoparticles) [sij . The presence 
of resonant transitions at lu and 3io in graphene not only 
enhances x^^\ but also ensm'es strong linear (and nonlin- 
ear) absorption at both frequencies. This absorption is a 
detriment to the third order processes, because it reduces 
the fundamental intensity needed for THG and leads to 
the subsequent absorption of the third harmonic [6^ . In 
fact, one can define a figure of merit (FOM) for the THG 
as the ratio of the total number of third harmonic pho- 
tons generated to the total number of photons (uj and 3a;) 
absorbed. It is straightforward to show that this FOM 
is, in turn, proportional to the ratio x^^V^^- Thus, while 
graphene and thin graphite films (and indeed all materi- 
als with resonant transitions) have large x^'^'''^: they do 
not have unusually large FOMs. For example, for our 
experimental conditions, a pulse with a fluence of 100 
/iJ/cm^ (peak intensity of ^ 0.3 GW/cm^) would have a 
third harmonic conversion efficiency of less than 10"^", 
but would lose 0.023 of its energy to linear absorption. 

In summary we have observed third harmonic genera- 
tion from graphene and few-layer graphite films. For a 
few layers, the emission scales with the cube of the in- 
tensity and the square of the number of atomic layers. It 
is characterized by an effective third-order susceptibility 
of the order of x^^^Sw; w, w) - lO^^^ m^/Y^. Such 
a strong third order response originates from the reso- 
nant nature of the light-graphene interaction, which also 



causes strong absorption of light. 
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